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SUMMARY 


Noise  data  for  normal  aircraft  operations  are  usually  cal¬ 
culated  by  the  NOISEMAP  computer  program  at  grid  points  1000 
feet  apart.  Data  from  blast  noise  and  supersonic  aircraft 
are  calculated  at  grid  points  several  thousand  feet  apart. 
This  report  describes  a  computer  program  that  was  written  to 
allow  the  two  sets  of  data  to  be  combined. 


PREFACE 


This  research  was  performed  for  the  Air  Force  Aerospace 
Medical  Research  Laboratory  at  Wright-Patterson  Air  Force  Base, 
Ohio  under  Project/Task  723107,  Technology  to  Define  and  Assess 
Environmental  Quality  of  Noise  From  Air  Force  Operations. 
Technical  monitor  for  this  effort  was  Mr.  Jerry  D.  Speakman  of 
the  Biodynamic  Environment  Branch,  Biodynamics  and  Bioengineering 
Division. 
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INTRODUCTION 

Noise  data  around  air  bases  are  normally  calculated  by  NOISEMAP 
at  1000  foot  grid  spacings.  However,  when  noise  from  supersonic 
aircraft  or  from  blast  data  is  evaluated,  it  is  often  desir¬ 
able  to  use  a  much  greater  grid  spacing.  The  purpose  of  this 
task  was  to  provide  a  means  of  adding  together  data  calculated 
at  different  spacings. 

A  general  technical  discussion  is  given  in  the  next  section. 
Operating  instructions  are  presented  in  the  final  section. 
Appendix  A  contains  a  listing  of  the  program.  Appendix  B 
contains  a  sample  run  including  a  plot  of  the  resultant  grid. 

TECHNICAL  DISCUSSION 

Each  set  of  noise  data  consists  of  10,000  data  points.  The 
distance  between  data  points  is  fixed  for  any  one  NOISEMAP 
computer  run,  but  may  be  varied  between  runs.  It  is  desirable 
to  calculate  the  noise  from  normal  aircraft  operations  using 
a  relatively  small  grid  spacing.  This  is  required  to  obtain 
a  sufficient  level  of  detail.  Noise  from  supersonic  aircraft 
or  from  blasts  may  cover  an  extremely  large  area  and  a  much 
greater  grid  spacing  may  be  used.  After  investigating  the 
two  separate  components  of  noise,  it  may  be  desirable  to  com¬ 
bine  the  data.  NOISEMAP  does  not  allow  different  size  grids 
to  be  added,  so  a  new  procedure  was  needed. 

After  investigating  several  methods  of  combining  the  grids,  a 
simple  approach  using  interpolation  was  found  to  be  the  most 
practical.  The  output  of  the  program  is  a  new  grid  of  10,000 
points  at  either  of  the  original  grid  spacings.  There  are  two 
basic  options  available  in  the  program;  the  small  grid  may  be 
added  to  the  large  grid  with  the  output  at  the  grid  spacing  of 

the  large  grid  or  the  large  grid  may  be  added  to  the  small  grid 

✓ 

with  the  output  at  the  grid  spacing  of  the  small  grid. 
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The  program  locates  the  small  grid  in  relationship  to  the  large 
grid.  If  the  two  grids  were  generated  with  consistent  coordinate 
no  modifications  are  required.  If  required,  the  origin  of  either 
or  both  of  the  original  dumps  may  be  redefined.  The  possibility 
of  allowing  one  grid  to  be  rotated  with  respect  to  the  other  was 
investigated.  The  coding  required  to  allow  this  would  be  much 
more  complicated  and  infrequently  used.  It  was  therefore  decided 
not  to  pursue  that  capability. 

The  program  first  reads  the  two  noise  dumps  produced  by  NOISEMAP. 
Each  dump  contains  a  header  that  identifies  the  following: 

.  Number  of  the  dump 

.  Logical  unit  on  which  the  dump  was  written 
.  Type  of  noise  data  (i.e.  DNL  or  NEF) 

.  Date  the  dump  was  written 
.  Title  from  the  • AIRFLD'  card 
.  Grid  spacing 
.  Field  altitude 
.  X  origin 
.  Y  origin 

.  Magnetic  declination. 

An  array  of  100  x  100  noise  data  points  follows  the  header. 
Finally,  additional  data  required  by  NOISEMAP  are  given. 

A  new  header  is  then  written  for  the  new  combined  tape.  Several 
items  may  be  changed.  The  dump  number  will  always  be  1.  The 
logical  unit  number  can  be  defined  or  will  be  14  by  default.  The 
date  will  be  the  current  date.  The  title  is  replaced.  The  grid 
spacing  will  be  one  of  the  two  original  grid  spacings  depending 
on  the  option  selected.  The  X  and  Y  origin  may  be  redefined. 
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If  the  small  grid  is  added  to  the  large  grid,  only  the  area  of 
overlap  is  affected.  If  the  location  of  a  large  grid  point 
coincides  with  the  small  grid  point,  a  simple  addition  is  re¬ 
quired.  If  the  points  do  not  coincide,  a  linear  interpolation 
is  performed  in  the  X  and/or  Y  direction.  Code  was  prepared  to 
extrapolate  the  data  outside  the  area  of  the  small  grid.  In 
testing  the  algorithms,  it  was  found  that,  at  times,  the  noise 
levels  were  increasing  at  the  boundary  and  would,  therefore, 
continue  to  increase  if  extrapolated.  The  purpose  of  the 
extrapolation  would  have  been  to  made  sure  that  the  transition 
appeared  smooth.  Instead,  an  area  of  uncertainty  was  developed 
so  the  code  was  removed. 

Adding  the  small  grid  to  the  large  uses  the  same  basic  routines. 
The  entire  small  grid  will  be  updated.  A  simple  linear  interpola¬ 
tion  is  used  between  large  grid  points  to  get  the  values  to  be 
added  to  the  small  grid. 

OPERATING  INSTRUCTIONS 

There  are  five  types  of  control  cards  that  may  be  required  to 
operate  the  program.  Each  type  has  a  keyword  and  is  discussed 
in  detail  in  the  following  pages.  They  follow  In  alphabetical 
order  but  may  be  inserted  into  the  input  deck  in  any  order.  Un¬ 
less  otherwise  defined,  it  is  assumed  that  the  large  grid  data 
are  stored  on  logical  unit  12,  the  small  grid  data  are  stored 
on  logical  unit  13  and  the  combined  data  will  be  written  on 
logical  unit  14. 

Of  the  five  cards,  only  two  are  required  for  the  program  to 
operate.  The  first  is  the  TITLE  card  which  provides  the  airfield 
title  to  be  written  on  the  header  to  the  dump.  The  second  re¬ 
quired  card  is  the  OPTION  card  that  determines  whether  the  small 
grid  is  added  to  the  large  grid  or  the  large  grid  is  added  to 
the  small  grid. 
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KEYWORD  BGRID 


7 . 5*5 . 5*3 . .  .7*5 


/bgrid 


12-  100000.  100000. 


? . iti5 . a 


3 . MB 1 


Function:  To  define  the  logical  unit  and  X  and  Y  origin 
for  large  grid  (optional). 


Columns 


1-6 

7-14 

15  -  22 

23  -  30 


BGRID 

Number  of  the  logical  unit 
containing  the  small  grid 
data  (optional).  Default 
value  is  12.  F  8.0  format. 

Redefinition  of  the  X  origin 
of  large  grid  (optional). 

F  8.0  format,  (feet) 

Redefinition  of  the  Y  origin 
of  large  grid  (optional). 

F  8.0  format,  (feet) 
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KEYWORD  CGRID 


CGRID 


If. 


7 . ins . i. 


s . 7*1  .  .7*$ 


Function:  To  define  the  logical  unit  for  combined  grid 
(optional) . 


Columns 


1-5  CGRID 

7-14  Number  of  the  logical  unit  to 

be  used  to  write  combined  grid 
(optional).  Default  value  is 
14.  F  8.0  format. 
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KEYWORD  FGRID 


KEYWORD  OPTION 


/  OPTION]  1. 


7 . 1*5 . 5! 


.  .5*5 


.7*1  .  .7*5  . 


Function:  To  identify  whether  the  small  grid  is  added 

to  large  grid  or  large  grid  is  added  to  small 
grid. 


Columns 

1-6  OPTION 


7  -  I*. 

Large 

grid  added 

to 

small 

grid . 

Enter  1. 

Small 

grid  added 

to 

large 

grid. 

Enter  2. 
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Function:  To  provide  the  alpha-numeric  data  to  be  used 
as  a  title  for  the  airfield. 

Columns 

1-5  TITLE 


7-76 


Alpha-numeric  data  to  be  used  for 
the  title. 


OUTPUT 


The  printed  output  from  the  computer  program  is  very  simple. 
It  identifies  the  tape  header  for  each  of  the  tapes  and 
identifies  the  option.  A  sample  of  the  output  follows.  A 
dump  of  noise  data  is  output  to  the  specified  logical  unit. 
This  dump  can  then  be  read  by  NOISEMAP  for  developing  printed 
grids  or  plots. 
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TAPE  HEADER  FOR  LARGE  GRID 


DUMP  1  UNIT  13  PROGRAM  DNL  DATE  03/29/7R 

AIRFIELD  EDWARD  AFB  SUPERSONIC  A/C 
GRID  SPACING  12000.  FIELD  ALTITUDE  2302. 

GRID  ORIGIN  X  -20000D.  GRID  ORIGIN  Y  -200000- 

MAGNETIC  DECL  If. 82 

TAPE  HEADER  FOR  SMALL  GRID 

DUMP  1  UNIT  13  PROGRAM  DNL  DATE  11/17/78 

AIRFIELD  EDWARD  AFB  SUPERSONIC  A/C  +  RANGE 
GRID  SPACING  2000.  FIELD  ALTITUDE  2302. 

GRID  ORIGIN  X  2920D0.  GRID  ORIGIN  Y  100000. 

MAGNETIC  DECL  If. 82 

TAPE  HEADER  FOR  COMBINED  GRID 

DUMP  1  UNIT  If  PROGRAM  DNL  DATE  0f/0f/7R 

AIRFIELD  EDWARDS  AFB  COMBINED  TAPE 
GRID  SPACING  2000.  FIELD  ALTITUDE  2302. 

GRID  ORIGIN  X  fR2000.  GRID  ORIGIN  Y  300000. 

MAGNETIC  DECL  If. 82 


LARGE  GRID  ADDED  TO  SMALL  GRID 


SAMPLE  OUTPUT 
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APPENDIX  A 


COMPUTER  PROGRAM  LISTINGS 


PSDbkAn  AOUokD  (INPUT,  j  UTPU  T .  TAP_S=InPUT»  TAPc 6  =  uUTPuT  , 
1  I  APE12*S1A,TAPE1  3-SIa.TAPE  IA  =  :>IA  > 

C 

C 

REAL  KfcYWO 

REAL  LAttEL 

LOGICAL  LbU.LFu 

LOGICAL  BGF.FGf 

EQUIVALENCE  (  I D  JN  P  ,  HE  AUt  r-.  (2  I  1 

ttiu  lVAi_tNCE  < lUNlT,H£AUbk( 3 >1 

C'JUNCN  /GRIDS/  N0F,  NBF 1.  »  isb  (  IGo  ,  100  I  ,  aS,  FGdOO.lUOl,  F 
1  bX»  JY.  kS*  rtL*  NLi  RMH.  HH,  k.Ni'1,  NN,  T 
CIMENS  IUN  I  TEXT!  7  >  ,  HE  ALe*  t  30  »  ,  LAbl.u  (  20  J  »KOh(  10 C» 

OintNSlUN  CARDIE) 

DATA  TITLE.  FjRIO,  DGRID,  OPTION,  CGklD 
1/  5H TITLE  .SHF GRID  .  SH8bklO,  bMuPTlUN,  SHCbklD  / 
DATA  LbG.LFL.FGF.3GF  /.FALSt.,. FAi.SE. ..FALSE*, .FAlSa./ 
DATA  T/2S.0 / 

DAI  A  Nl/5/.  NC/O / 

DATA  NaG.NFG.NLG  /  12.  13.  1A  / 

DATA  NbF  /100/ 

C 

L 

C  IFLAG  IS  US tD  TU  INuICaTL  IF  USck  » ISHtS  TD  GG  FKGn 

C  aiG  v»k  I D  TO  FINE  bRlD  (IFLAG*lI  JK  FKGn  FINE  bRlD  TO 

C  cIG  bK ID  (  IFLAG *2  >. 

C 

KEklNO  NBC 
REnIND  NFG 
CALL  GaTEUdATLI 

1  kEAO(S.SOOU)  C AkO 
IF(SGFtSl)  9999,2 

2  CON  TINGE 

JECUl>£(i0»5103,CARDt  KEYwG 
IF  (  KEYWD  -  TITLE!  3,3b 
j  IF  <  KEYWD  -  FGKIO  »  A, AO 

A  IF  I  KtYWQ  -  BbkIO  »  b,»U 

5  IF  (  KcYWO  -  OPT  I  ON  I  6,60 

6  IF  (  KtYWO  -  CGkIO  I  7,/0 

7  WP1TF  (  6.60111  KEYWO 
STOP  l 

30  DECODE  170,6101, CARD!  IlASEl t I » , I =b , IS  I 
GO  TU  1 

AO  DEL  DUE  t  70,  S  102  ,  CAR  0  )  XFG.FXO.FYO 
IF  (XFb  .EQ.  0.)  GU  TO  AL 
NFG*XFb 
AL  CGnTINUc 

ififxo.ne.o.i  gd  TU  as 

TEST  *  3  I GN  < 1 • ,F  XO ♦ 

IFlIEST)  l.AS.AS 
AS  lFd  «  .TRUE. 

DO  TU  i 

SO  DECODE  (  70, 6102, CARD)  XuG.bXO.jYD 
I  F  (  XBb  .EU.  0.1  GU  TO  V2 
Nbb»XBb 
Sc-  CONTlNUt 

If  ( aXO.NE .0. )  Go  TD  SS 
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TfcjT  -  6IGN(  i.,dXOI 
IF(TESI)  1  ,66.66 
66  LdU  -  .TRUE. 

CO  TO  1 

60  OtCUOE  (  70. 6102.CARU)  Xr  LAD 
I FlAG-XFLAG 

GO  TO  1 

70  DECODE  (  70.6  1 02. CARD  )  XCG 
NCG-XCG 
GO  TO  1 
9999  CONTINUE 

I F ( I  FLAG. £0. 1 )  FGF-.TRUt. 

IF! 1FLAG.EQ.2)  8GF  *  .TKot. 

100  REAO(NBG)  HEADER 

IF (EOF (N8G  )  I  160.110 
110  I F ( He ADE R(1).EU.3HEND)  uU  TO  16. 

PRINT  6000 

PRINT  6010.  (HEADER! 1 1 ,i -1,201 
8S-  HEADER! 16  ) 

I F (LBO  I  GO  TO  III 
bXO-  HEADER (101 
B  YO  »  HE  AOER ( 19) 

1 1^  CONTINUE 

HEADER! 5 )  -  TD AT t 
HEADER ( 18  )  *  BXO 
HEADER! 19}  *  6Y0 
IDUHP-1 
IUNIT-NCG 
00  11J  1-6.16 
113  HEADER ( I  I  »  LABEL!  I  I 

IF  ! BGF )  WR1TEINCG)  HEADER 
IF(BGF)  PRINT  6002 

IF ! BGF I  PRINT  6010. (HEADER! I ) .1-l.cOI 
DU  120  J-1,100 
REAO(NBG)  !6G(  1.  J  >  »  I  - 1.  j.  001 
120  CONTINUE 
160  READ (NFG  >  HEAOEk 

I F ( EOF (NFG I  )  200.160 
160  IF (HEAOE R ( 1 1  .EG.  3HENJ)  GJ  TO 
PRINT  6001 

PRINT  6010  .  (HEADER! II .1-1.20) 

FS  -  HEADER ( 16 ) 

IF(LFO)  GO  TO  162 
FXO  -  HEAOER(ld) 

FYO  -  HEA0ERO9) 

162  CONTINUE 

HEADER ( 6  I  -  TDATE 
HEA OER ( 1 8 1  -  FXO 
HEADER! 19)  -  FYO 
I  UNIT -NCG 
ICUNP-1 
00  163  1-6.16 

163  HEAOER! I  I  -  LABEL (I  I 
IF1FGF InRITE! NCG I  HEAOER 
IF (FGF I  PRINT  6002 

IFIFGFI  PRINT  60i0.(HEAuER( 1 1 . 1-1.20) 
00  170  J-1,100 


16 


170 

200 


C 

c 

0 


300 


350 


360 


370 


0 

c 

c 

r 

v 

500 


AGO 


A  1 0 


A20 


■>000 

5101 

5102 

>>103 

6000 

6001 

6002 

o003 

oOOA 


KtAolMFG)  (  F  C  ( 1  •  J  )  » 1  *  1  **.  00  I 

CUNT  1 Nut 

CONTINUE 

call  calc  c b xo»  aro,  fxo»  fyo> 

IF  (IFLAC.E0.1I  GO  TO  5oO 


USc  FINE  GklO  TU  GET  Blu  GRID  VhLul  S 


PRInT  6005 
PRINT  oOQA 
PRINT  6005 

IF  (  (MM.EQ.MLI  .ANJ.  (NN.EO.NLU  o..  TO  jOO 

call  inner 

CONTINUE 

u'J  350  J-1,100 

RRiTE(NCG)  ( BG( I  * J)  *  1*1, 10J  I 

C On  T 1  NOE 

S.cAUINLlI  ITEXT 

wKITECNCCI  (TEXT 

NAnn  =  iT£XT(l» 

IftNANNI  36U,5G 
O'Z  370  i  *  1 »  NANN 
ncmlMNbGI  ITEXT 
hkiTEINCG)  1  TEXT 


Cun T InUl 
ST  up  2 


uSc  dlu  GRIO  TO  GET  FIn:.  GRID 

CGNT INUE 
print  6005 
PRINT  t»OC 3 
PRINT  o005 
C Ai_L  dOTFG 
UU  AOO  J*1,10C 

WKITE(NCG)  ( FG( i «  Ji  *1*1,1001 

CONTINUE 

READ(NFGl  I  TEXT 

NR  1 T  £( NCGI  ITEXT 

N  Ann  =  ITExT ( 1  ) 

IF  (NANN)  A  10, 50 
Du  A 20  1 =1 »  NANN 
KtAU(NFG)  ITEXT 
RKi Tfc ( NCG  I  ITEXT 
CUNT InuE 


VALUi  5 


& 


& 

•N' 


<-v 


<6^  ov 
* 


STOP  3 


FORMAT  ( 8 A  10  ) 

FORMAT  ( 6X i 10 Ac ,  AX  I 
FORMAT  (6X,aF8.0) 

FOnmA  T( A6« A  X  > 

FORMAT  (  2  7iliT  APt  HEADER  FUr  LA  Kut  „*I0  ) 

FORMAT  ( 2  7HGTAPE  HEADER  F0«  SMAlL  *Klu  ) 
FORMAT  (  30HG  TAPE  HEADER  FOR  CUMoiNl.U  GnI  0 
FORMAT  (31M0LARUE  GRID  ADDED  TO  SH„lL  GkIu 
FORMAT  (  3 lHOSHAbL  GklD  m ODr.  0  Tu  L  A*.  uE  GRID 
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6005  FOiirtATUlHO********  ♦♦****♦*•*♦♦•**♦•*•*♦•**«***♦**♦  » 

6010  F  OKMAT  ( 1H0.  10X»  A*  *  I*.  6H  UNIT  13.  9H  PROGRAM  A  7,  6H  UAlt  UO/ 

l  1H  »  10X»  9HAIRF IELD  10A6  / 

i  1H  ,  10X,  13HGR10  SPACING  h9,0.  bX.  IbrtFltLO  ALTITUDE 

H  F7.0  / 1H  •  10X.  1AHGRIU  ORIGIN  X  Ftt.O,  5X, 

5  1AHGR1D  ORIGIN  Y  Fb.O  /  Ih  .  10a.  UHMAGNETiC  Gt  Cl  f-b.i  I 

6011  FORMAT  I 19H  INVALID  KEYWORD  *A1GI 

END 


SUokuUT  1  NE  CALC  <  BX  0  *  b  YO»  FXO.  fYi,  I 

COMMON  /GRIOS/  NBF  »  N0F1 »  bG (  IGu  .  lu  U  )  *  bS.  H. t 100  * 1.00  I  *  FS * 
i  OX.  uY*  RS.  ML.  NL,  Rum,  MM,  kNNi  nn.  T 
c 
c 

C  NBF  *  NUMBER  OF  ROWS  ICOlUMNS)  iN  frtE  UR  I OS 

C  NBF 1  =  NBF  -  1 

C  BG  *  AkRAY  UF  BIG  GRID  POINTS 

C  BS  *  SPACING  IN  BIG  GRIu 

C  d  XU  *  b Yu  -  X  »  Y  CuOPD I NATc S  FUR  UrIGIn  OF  bIG  Gk  lu 

C  FG  -  ARRAY  UF  FINE  GRIO  VALUES 

C  F  S  *  SPACING  FUR  FINE  GkIO 

C  FXO.FYO  -  X  » Y  CUURDINATcS  FOR  U*lGlN  OF  Fl.sfc  GRID 

C  RS  *  BS  /  FS 

C  OX  -  FXO  -  8X0  ♦  I DX-1 ) ♦ ; } 

C  OY  -  FYO  -  BYO  ♦  ll)Y-l»*BS 

C  ML.NL  -  I  NO  ICE  S  OF  LOWEST  X.Y  COORDINATES  UF  ThOSE  BIG  GKIc 

c  points  contained  in  fi Nt  grid 

C  MM . NN  -  INDICES  OF  HIGHcST  a.Y  COORDINATES  bF  ThUSt  BIG  Gki^ 

c  points  contained  in  fine  grid 

c 

c 


NBF  1 

*  NBF  —  1 

LX  * 

(  FXO— 8X0 1  / 

BS  ♦  1. 

DY  = 

(FYO— 8Y0I  / 

as  ♦  i. 

RS  * 

BS/FS 

ML  * 

IFIXIUX) 

NL  * 

IF  IXIDYI 

R9-#  *  N8F1  /  RS 
KMn  «  uX  ♦  K99 
R  NN  *  OY  ♦  K99 
MM  =  I F I X  I R  MM ) 

NN  -  IFIX(RNN) 

IF  ( DX— ML  • NE •  01  ML  -  ML  ♦  I 
IF  IDY-NL  .NE.  U.l  NL  «  NL  ♦  l 
RETURN 
END 
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o  nr 


SUBROUTINE  INNER 


C 


C 


*00 

300 


INNER  luops  on  all  BIG  GRIo  POINTS  CUNTAINLO  in  THE  FlNt  Gn i 
INTERPOLATING  FINE  GRID  POINTS  Tl  GET  VALUES  FOR  lACh. 

CUrtNON  /GRIGS/  N8F,  N8F1,  oO { lOu , 1 uO »  ,  oS,  FGIi0C,100l,  FS, 
l  GX,  OYf  RS,  ML,  NL»  Rrltt ,  rtM,  *NN,  NN»  T 

yj  *  (nl-oy-1.1  *  r s  ♦  i. 

00  t>00  J*NL  * NN 
YJ  *  YJ  ♦  RS 
JV  *  1FIXIYJI 
A I  *  I HL-UX-1. I  *  RS  ♦  1. 
ju  *00  I»HL*«rt 
XI  -  XI  ♦  RS 
IX  *  I F I  XI X I ) 

CALL  INTERFI IX, JY,X1 ,YJ  ,GI 
6G  < I , J )  *  UO  < I , J  I  ♦  G 
CONTINUE 
CONTINUE 
RETURN 
ENO 
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n  r>  0000000  r: 


SUBROUTINE  INTERF  (l*  J*  RI  ♦  RJ*  Ul 


INTERF  INTERPOLATES  F iNt  GRID  POINTS  TU  GET  A  SlNGLt  VALUE 
FUR  A  BIG  GRID  POINT. 

I*J  ARE  THE  1N0ICES  FOR  THE  FINE  GRlU  POINT  CLOSEST* 

BUT  TO  THE  LEFT  ANO  6ELGW*  The  UcSlREU  blG  GRID  POINT. 
KI.RJ  ARE  THE  ACTUAL  FLOATING  POINT  CuOROINATES  THE  BIG 
GRID  POINT  WOULD  HAVE  WERE  IT  IN  THE  FINE  GKIL*. 

COMMON  /GRIDS/  N BF  »  NBFi *  uG ( 100 , 10 0 1 ,  BS*  FG(100*100I*  FS* 
i  OX*  OY*  RS*  ML*  NL*  RMM*  MM*  KNN*  NN*  T 

INT ER PUL ATE  IN  X  DIRECTION 

IF  < R I  .EO.  I)  GO  TO  20 0 

GX  *  (PGI  I  ♦  I  *  J )  -  FGtl»J>>  *  tki-ll  ♦  FuU.Jl 
if  (RJ  .NE.  Jl  GO  TU  *00 
G  *  CX 
RETURN 
C 

C  NO  INTERPOLATION  NiiEDEO 

C 

200  CONTINUE 

IF  (RJ  .NE.  J»  GO  TO  30u 
G  =  F  G ( I  »  Jl 
return 
c 

G  INTERPOLATE  IN  Y  OIkECTIUN 

0 

3o  0  CONTINUE 

G  *  (FG(I*J*1I  -  FG(I*J)I  *  (RJ-JI  ♦  FG(1.J» 

RETURN 

c 

C  INTERPOLATE  IN  oOTH  DIRECTIONS 

C 

*  u  0  C  ONT  I  NGc 

GY  «  (FG(I+1«J+1I  -  FG( I , J* 1 II  *  (vi-II  ♦  FG(I*J*L> 

G  =  (GY  -  GXI  ♦  (RJ-JI  *  GX 

Rt  T  GKN 

lng 
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m  nn  r-  n  n  f)  onr»nrtf>nnnnn 


SUBROUTINE  INTERS  (  V L  ,Vc  ,D.  R  ,0E  ,  G ,  ~i  G  > 

INTERB  INTERPOLATES  BIG  GRID  PolNTS  Tli  GET  VALUES  FUR 
FINE  GRID  POINTS 

VI,  V2  -  GRID  VALUES  TU  BE  INI; RPOlA TEG 

D  -  DISTANCE  (IN  BIG  GRID  SPACING)  TU  FIRST  FINE  GRID  POINT 
FROH  VI 
R  -  FS  /  BS 

OE  -  OUTPUT  VALUE  JlSTANCE  FROM  V2  TO  NEXT  FINE  GRID  POINT 

G  -  ARRAY  OF  INTERPOLATED  VALUES 

NG  -  NUMBER  OF  VALUES  IN  G 

DIMENSION  G(l) 


C  *  V2  -  Vi 
G  ( i  I  *  VI  ♦  C*C 
C  *  C  *R 

IF  (NG  .NE.  01  GO  TU  30u 
NG  *  1 
DE  *  D 

IF  (OE  .GE.  II  RETUKN 

0  CONTINUE 

tit  *  OE  ♦  R 

IF  IDE  .GT.  i. »  Rt  TURN 
NG  *  NG  ♦  1 
G (NGI  -  GING-II  ♦  C 

GU  TG  100 


00  CONTINUE 

IF  (NG  .EG.  1)  RETURN 
DO  *00  1*2, NG 

G(  II  *  G( I  —  1 1  ♦  C 
*00  CONTINUE 
RETURN 
END 
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n  o  n  non  n  o  c*  ci 


SUBROUTINE  OGTFG 


C 


oGTFG  USES  BIG  GRIG  PGanTS  TU  INTERPOLATE  VALUtS  PUR  *uL  Hwt 
GRID  POINTS. 

CONHON  /GRIDS/  NBF,  NBFl»  jG  4  IO'j  .  ID  0  )  *  cS.  P  G(  100  *  1  u  0  )  ,  PS. 

1  DX.  Ur,  k  S ,  ilC.  NL,  RrtN,  (IN.  *NN,  NN ,  T 
COrtMON  /SCRACH/  GK100I.  G<»  (  LOG  I  .  ,j(lOOI 

INITIAL  COMPUTATIONS 


NNL 

=  ML  -  1 

IF 

( ML  .EQ.  DX) 

NHL  * 

Hl 

NNL 

n 

2 

r 

t 

*— 

IF 

(NL  .EQ.  DY) 

NNL  * 

Nt. 

11  MM 

=  .'IN  ♦  1 

IF 

(NH  .EQ.  RNN) 

NNN 

*  in 

NNN 

*  NN 

IF 

(NN  .EU.  RNN) 

NNN 

*  N  N 

HP  1 

=  NHL  ♦  1 

XOfc 

*  0. 

voe 

*  0. 

XOO 

*  OX  -  NHL 

ro 

*  CY  -  NNL 

KRS 

«=  l.  /  RS 

L  JO 

*  0 

LOOP  ON  BIG  GRID  POINTS  THAT  SUakUjNO  PINE  GkIl 


DO  300  J-NNL.NNN 
NO  L  *  0 

GALL  I NTE*B(8G(MML» JI.3GI H«L» i* i) . YO *RRS • YGfc t G 1 »NG l ) 

NG2  *  NBF  -  L JO 

NG1  ®  M I  NO ( NG1 . NO  3  I 

NG<2  «  NG1 

XD  *  XOO 

LIO  »  0 

it  U  600  I  *M  P l .  HMM 

CALL  INTERU  (6G(I.JI.8G(  I  »  J*  1 1  ,  VO  .Rk  S  *  YDF  .LX  »NG/ » 

N  «  NBF  -  LiO 
NC3  «  0 
LJ  »  LiO 
DO  <*00  K*  1»NG<? 

CALL  I  NTEK8  (  G  l(Kl  .Gl(K).Au.kkS.XlE.G3.NGJ) 

LJ  *  LJ  ♦  1 
NG 3  «  N I«0 ( N  »NG  3 1 
LI  -  LIO 
DU  300  L-l.NGJ 
LI  *  LI  ♦  I 

FG(Ll.LJ)  «  FG(ul.LJ)  *  G  j( w  I 
300  CONTINUE 

<*00  CONTINUE 

OU  *j 00  X*  1.NG2 
G1(KI  «  G3 ( K I 
SOO  CONTINUE 
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SUtrtCUTJNE  tiGTFo 


?<*/  7*  UF  T  *i 


oO 


t)5 


xo  *  xot 

LiO  *  LI 


60  0 

tONTINUt 

4  JO  -  LJ 

*i>  *  YOE 

doo 

COn I  I Nu£ 

kEIuKn 

t  NL 
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-  IMXIXDtl 

ifixir jei 


APPENDIX  B 


SAMPLE  COMPUTER  RUN 
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TITLE  EDWARDS  AFB  COMBINED  TAPE 
OPTION2- 

BGRID  12.  0.  0. 

FGRID  13.  *<12000.  3D0D00. 


EXAMPLE  INPUT  DATA 


TAPt  HEADtR  FOR  LARGi.  CR  it) 


DuMP  l  UNIT  13  PROGRAM  UNL  DATE  0J/29//9 

AIRFUlO  EUWARu  AFb  SUPERSONIC  A/C 
GRID  SPACING  12000.  FIELD  AlTITuGE  23u2. 

GRID  ORIGIN  X  -ZOOuOO.  GRlb  ORIGIN  »  -^OGO^O. 

MAGNETIC  DLCL  14.62 

tape  header  for  combined  grid 

DUMP  1  UNIT  14  PROGRAM 
AIRFIELD  EDWARDS  AFB  COMBI 
GRIO  SPACING  12000. 

GkIU  ORIGIN  X  0. 

MAGNETIC  DtCL  14.62 

T APt  HEADER  FOR  SMALL  GRID 

DUMP  1  UNIT  13  PROGRAM  UNL  u ATE  li/17/7e 

AIRFIELD  EDWARD  Af  B  SUPERSONIC  A/C  ♦  RANGE 
GRIO  SPACING  2000.  FIElO  AlTITuGE  23u2. 

GRIO  ORIGIN  X  292000.  GRID  ORIGIN  T  IOOOlO. 

MAGNETIC  Ot CL  14.62 

********************* ******************* 

SMALL  GRIO  ADOEO  TO  LARGE  oRlD 

********r******************************* 


UNL  DATE  04/lB/oU 

NED  TAPE 

FIELD  ALTITUDE  2 3» 2 • 
GRID  ORIGIN  Y  u. 


EXAMPLE  OUTPUT  DATA 


SAMPLE  COMPUT 


